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THE EQUIVALENCE OF THE CONDUCTION THEORIES OF 
RASHEVSKY AND RUSHTON 


ALVIN M. WEINBERG 
THE UNIVERSITY OF CHICAGO 


The one-factor, semi-formal nerve conduction theory of Rash- 
evsky and the dynamical theory of Rushton are shown to be com- 
pletely equivalent. 


In 1933 N. Rashevsky proposed a semi-formal theory of nerve 
conduction based on the excitation theory of H. A. Blair (1932) and 
on the theory of local re-excitation by action currents. This theory is 
“semi-formal” in the sense that in it the excitation process is consid- 
ered purely formally, while the distribution of currents in nerve is 
calculated on the dynamical assumption that a nerve fiber is a pure 
resistance network (Figure 1 a). 

Four years later, W. A. H. Rushton (1937) proposed a dynamical 
theory of excitation and conduction in which a nerve fiber is consid- 
ered to be a resistance-capacitance network (Figure 1 b) and excita- 
tion occurs when 6, the charge (or voltage) across the nerve mem- 
_ brane (or capacitance in the network), exceeds a certain critical val- 
ue, 6. Excitation in both theories is characterized by the appearance 
of an action E.M.F. in series with the nerve sheath (Figures 1c, 1 d) 
at the excited point. 

We shall now show that the two theories are formally equivalent 
and that Rushton’s model may be considered as a dynamical instance 
of Rashevsky’s semi-formal theory. 

The fundamental equation of Rushton’s theory describes the vari- 
ation of charge, 9, along a line consisting of an unexcited portion 
(Figure 1 b) joined to an excited portion (Figure 1d). This equa- 
tion is given by Rushton (1937) as 


7 iN ROT RRL PRN Se) RR od 
Ge aa) Rm Raa, ca 
where 
r, = resistance per unit length of interstitial fluid. 
r, = resistance per unit length of axis cylinder. 
= = conductance per unit length of sheath for radial currents. 
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1 ; 1 1 
— = conductance per unit length of membrane. Ce eyes R? ; 


p 
6 = charge per unit length of membrane. 


¢ = capacitance of membrane per unit length. 
| dl Sy A a= Ue 
a = Cp 
I = i, +i, =sum of currents flowing parallel to nerve outside 
and inside. (I = 0 except between stimulating elec- 
trodes. ) 
E, = 0 in resting region. 
= F# in excited region. 


The derivation of Rashevsky’s fundamental equation requires two 
steps. First, an elementary application of Kirchhoff’s laws to a line 
consisting of unexcited and excited portions (Figures la and Ic) 
joined at the point of excitation gives for the space variation of the 
current, 7, flowing through the myelin sheath, R, 


2 2 
ot w Lay onttaten tle Undies Be ; (2) 


where the constants have the same meaning as before. Second, the 
time variation of the formal excitatory state, p, is given by Blair’s 
equation 


where K and k are constants, and excitation occurs when p exceeds a 
critical value p.. Elimination of 7 between equations (2) and (3) gives 


Lae 1 r, ol old Os 
cD Gt) =-$ Se (4) 

Equation (4) is identical with (1) if p/K is identified with the 
membrane charge 6, and 1/k is put equal to a, the relaxation time of 
the membrane condenser. 

Since all the results of Rushton’s theory follow directly from 
equation (1), and since Rashevsky’s theory leads to a formally iden- 
tical equation, it may be concluded that the two theories are identical. 
For example, the original expression for the velocity of the impulse 
was given by Rashevsky as 


= LIK (5) 
where 4, the safety factor, is defined by 
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he action current — rheobasic current 
rheobasic current : 


or, since current and voltage are proportional, 


__ action voltage — rheobasic voltage 
4 
rheobasic voltage 


Rushton’s velocity formula is 


ee (6) 
a 


which of course is identical with (5) if k = 1/a. 

Whether Rushton’s resistance-capacitance model is the only pure- 
ly electrical model which is equivalent to Rashevsky’s semi-formal 
model cannot be ascertained a priori. The development of other dy- 
namical models which are equivalent to these two may be of interest 
because they would allow for a greater flexibility in the interpreta- 
tion of conduction phenomena. 

This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbot Memorial Fund of the University of Chicago. 


a. b. 
RASHEVSKY MODEL CRESTING STATE) RUSHTON MODEL CRESTING STATE) 


c 


RASHEVSKY MODEL (EXCITED STATE) RUSHTON MODEL CEXCITED STATE) 


FIGURE 1 
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FURTHER CONTRIBUTIONS TO THE THEORY OF CELL 
POLARITY AND SELF-REGULATION 


N. RASHEVSKY 
THE UNIVERSITY OF CHICAGO 


In continuation of a previous study a case of self-regulating 
polarity is investigated, in which a metabolite neither enters nor 
leaves the cell, being produced in one part of it and consumed in an- 
other. Such a cell possesses a larger degree of autonomy and inde- 
pendence of external conditions than the previously discussed one. 
Application to cellular locomotion is indicated. 


This study is a continuation of a previous paper (Rashevsky, 
1940; hereinafter referred to as loc. cit.). 

Consider a spherical cell of radius ” , inside of which a substance 
is produced at a rate of gq gm cm™ sec”. Let the same substance break 
down and therefore disappear at a rate proportional to its concentra- 
tion c, and given by — be gm cm* sec”. 

Let the membrane of the cell be completely impermeable to that 
substance, so that the latter neither enters nor leaves the cell. Under 
these conditions a stationary state is approached, in which the rate 
of production q of the substance is equal to its rate of consumption 
be. This gives c = q/b. The concentration is uniform throughout 
the cell. 

Consider now, as in the previous paper (loc. cit.) that the cell 
contains colloidal particles, which act as an inhibitory catalyst upon 
the rate of formation q of the substance. Denoting by 1 the concen- 
tration of the colloidal catalyst, we shall again use equation (24) of 
loc. cit. For a constant c, n is also constant throughout the cell. 

Let now, due to any disturbance, a small increase in 1 arise in 

some region of the cell, producing thus an asymmetry. In this region 
q will decrease and therefore ¢ will decrease. This will result in a 
negative gradient of ¢ in the direction of that region and therefore in 
an accumulation of colloidal particles in that region with a resulting 
further decrease of c. As in loc. cit. we thus will find an asymmetry 
or polarity automatically established and maintained. The conditions 
for stability are the same as before. 

The calculations proceed very much as in loc. cit. Denoting by 
the subscript 1 the quantities referring to one hemisphere and by the 
subscript 2 those referring to the other hemisphere and remembering 
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that in the stationary state the total amounts produced in one hemi- 
sphere equal the total amount consumed, plus the total amount dif- 
fusing into the other hemisphere, we have 


Zn °Q1 — 311° bC, — ar,D (¢,— C2) = 0, (1) 


27 9°3o — Bnto*bC2 + aroD (¢,— C2) = 0. 
Solving for ¢, and c, and calculating the quantity ¢, — ¢, we find: 


Thee 


7.26 + 3D (@1 — Qe) - (2) 


Gq —G = 
For the relation between the quantity x = c, — ¢. and 2% — ”, we 
have equation (23) of loc. cit. In its dependence on gq, — qs equation 
(2) is similar to equation (18) of loc. cit. Hence, putting 
=A (3) 
we shall have for 7 = c, — ¢, formally the same expressions as those 
obtained previously. 

In view of the assumed impermeability of the cell membrane to 
the substance in question, such a cell possesses an autonomy and inde- 
pendence of external conditions to a much greater degree than before. 

An interesting consequence may be considered here. Let the mem- 
brane, while impermeable to the metabolite, be freely permeable to 
water. The osmotic pressure due to the concentration of the meta- 
bolite and directed outward, will be greater at the pole, at which the 
concentration is greater. Therefore a net resultant force in the di- 
rection of that pole will be present, and the cell in that case will move 
forward with the pole of higher concentration. At first sight this con- 
clusion sounds like a paradox, for we seem to have a movement due to 
internal forces. This however is not so, since inside the cell the sub- 
stance flowing from one pole to another exerts a diffusion drag 
(Young, 1939) upon the water, which even inside the cell remains 
stationary, due to the assumed complete permeability of the cell mem- 
brane to water. 

Somewhat similar cases of spontaneous locomotion have been 
studied by Gale Young (unpublished) for the case of an asymmetric 
flow of a substance from the cell. In that case the mechanism of lo- 
comotion is much like that of a rocket. In the present case the me- 
chanism may be to some extent compared with a sort of caterpillar 
drive. Such mechanisms may well be responsible in some cases for 
spontaneous locomotion of unicellular organisms. 


x 


N. RASHEVSKY 67 


The author is indebted to Mr. H. D. Landahl for checking the 
algebra. 

This work has been aided in part by a grant from Dr. Wallace C. 
and Clara A. Abbot Memorial Fund of the University of Chicago. 
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CONTRIBUTIONS TO THE MATHEMATICAL THEORY 
OF ORGANIC FORM. II 


ASYMMETRIC METABOLISM OF CELLULAR AGGREGATES 


N. RASHEVSKY 
THE UNIVERSITY OF CHICAGO 


The case is studied of a cellular aggregate forming a hollow 
spherical shell, the physical constants and the metabolic rate being 
non-uniform throughout the shell. This leads to asymmetries of con- 
centration distributions, which are calculated here for some simple 
cases. 


In a previous paper (Rashevsky, 1940 b) we have investigated 
the mechanism of formation of a hollow shell from a densely packed 
cellular aggregate. This has been done especially with a view of pos- 
sible applications to the formation of a blastula. 

When, as is actually the case, the original cell from which all the 
cells of the blastula descend, is not uniform in its physico chemical 
constitution, then the cellular shell, forming the blastula, will also be 
nonuniform. This nonuniformity may manifest itself in a variation 
of the diffusion coefficients and permeabilities along the surface of 
the shell or in similar variations of the rates of reactions or both. 
For reference in future studies, we shall derive here some formulae 
for such 2 nonuniform shell. For simplicity we shall consider here 
the case of an infinite surface permeability and of an infinite external 
diffusion coefficient for one substance, produced or consumed by cells 
of the shell. 

Let the approximately spherical shell have an internal radius, 
7%), a thickness 6 < < 7 and be characterized by the following non- 
uniformity. Within a gingle-connected region 1, occupying a fraction 
a of the surface of the shell, the diffusion coefficient will be assumed 
to be D,, the rate of reaction q,. In the remaining part of the shell, 
2, the diffusion coefficient has the value D, , the rate of reaction — the 
value g.. The external concentration of the metabolized substance is 
denoted by ¢. Due to the infinite permeability, ¢) is also the inside 
concentration of the substance at the outer surface of the shell. 

If the shell were uniform in all respects, then the concentration 
inside the shell at the inner surface would also be the same at all 


BO an 
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points of the surface and within the cavity. Due to the nonuniform- 
ities assumed above, we shall have in general different concentra- 
tions ¢, and c in the two regions. Again, due to the assumption of 
infinite permeability, c, and c. represent also the concentrations of 
the substance inside the cavity in the immediate neighborhood of the 
shell. 

Due to this nonuniformity of concentration inside of the cavity, 
there will be not only a transport of the substance between the shell 
and the external medium, but there will also be a flow of substance 
inside the cavity. This flow will be determined by the average gradi- 
ent and the cross-section through which the flow takes place. The 
average gradient in the direction from 2 to 1 is of the order of mag- 
nitude of (c, — c.)/d where d is the average distance over which ¢ 
varies from ¢c, to c.. Therefore, denoting by D, the diffusion coeffi- 
cient of the liquid in the cavity, the flow per cm? per second is of the 
order of D,(c,; — ¢.)/d. The assumption of a finite D, and infinite 
external diffusion coefficient D, is not a self-contradicting one, for the 
cavity may contain solutes, for which the cells of the shell are imper- 
meable and which are practically absent in the outside. In any case 
the assumption D, = o does not introduce any fundamental restric- 
tions. The total flow within the cavity from region 1 to 2 is equal to 
D,(¢, — ¢2)S/d , where S is the average cross-section, through which 
the flow takes place. 

From physical considerations it is apparent that when a is very 
small, then practically everywhere inside the cavity the concentrations 
ce will be equal to c, , and only in a very small region will ¢ be equal 
to c,. Therefore the average gradient will extend only over a very 
small region. Thus with a tending to zero, d also tends to zero. It 
also tends to zero as a tends to 1. For a = 4, when just half the 
shell has one set of constants, and half another one, d has a maximum 
and is of the order of magnitude of 7. Approximately we may there- 
fore put 


d = 4r,a(1—a) , (1) 


which is one of the simplest expressions having the required property. 
Similarly the average cross-sectional area S is zero for a = 0 
and a = 1, and is of the order of ar,? fora = 4. We put 


S = 4nr,2a(1—a). (2) 
With these approximations the flow within the cavity is equal to 
AY Dy (c, a Co) e (3) 


The concentrations c, and c, are determined in the usual] way by con- 
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siderations of material balance (Rashevsky, 1940 a, chapter i). The 
amount of substance metabolized in the region 1 equals 4zr,? adq,. It 
is also equal to the difference between the amount 4zr,? aD, (¢: — €o) /6 
flowing outward and the amount given by (3), coming to the region 1 
from within the cavity: 


i 
Anr,? oD, an A-nt)Do (C1 — C2) = Ants? 050; « (4) 
Similarly, for region 2: 


Aree a), ae nrg (0p 6s) dare (1: —a) 602. 
(5) 


r) 
The second term, on the left hand sides of equations (4) and (5) have 
opposite signs, because when inside the cavity the substance is trans- 
ported towards region 1, it is transported away from region 2, and 
vice versa. 
Introducing the notations 


D=oaD,+ (1—a)Dz ; q=oq.+ (1—a)@, (6) 
we obtain from equations (4) and (5) after some rearrangements 
aos 2 37 
C= ae 4a(1 a) Dd Toi a Pod qd (7) 
4a(1—a)1%D,D, + DoD 6 
and 
4a (1 — a) D,6?7oq2 + Dod°q 


= (8) 
4a (1 —a)%DiD.+ DD 6 


Co = Cot 


For D,=D.=D; a= =q equations (7) and (8) reduce to 


6q 
C, = CG = C+ 7 y (9) 


an expression that can be obtained directly. For not too small values 
of a and for D, ~ D, ~ D, equations (7) and (8) reduce to 


6°q1 0°q2 
= ce a ; 1 
Cpesgee Do Gee Tame (10) 
The author is indebted to Mr. H. D. Landahl for checking the 
algebra, and for valuable suggestions and comments. 
This work has been aided in part by a grant from the Dr. Wallace 
C. and Clara A. Abbott Memorial Fund of the University of Chicago. 
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A CONTRIBUTION TO THE MATHEMATICAL BIOPHYSICS OF 
PSYCHOPHYSICAL DISCRIMINATION III 


H. D. LANDAHL 
THE UNIVERSITY OF CHICAGO 


i The neural mechanism previously discussed is further general- 
ized. The case is considered in which a random variation is asso- 
ciated with each stimulus. The mechanism is generalized and equa- 
tions are derived for discriminations between stimuli differing in 
several modalities. The latter indicates an analysis by the factor 
method. Suggestions are made in connection with the use of triads 
and with the problem of a multidimensional psychophysics. 


It is the purpose here to generalize, in two different directions, 
the neural mechanism previously described (Landahl, 1938, 1939; 
hereinafter referred to as A and B respectively). Reference shall be 
made to the diagramatic representation of the mechanism of Figure 
1 of the first paper (A) by writing Figure 1 (A). The notation used 
here shall conform to that previously used and is largely the same 
as the notation used by N. Rashevsky (1938). 

I. One modality but with independent fluctuations. In the pre- 
vious discussions it was considered that a random fluctuation occurred 
only at one synapse since the results would be essentially unchanged 
by introducing independent fluctuations at both s, and s. of Figure 
1(A). We proceed to show that this is so, then to investigate some 
conditions under which the more general consideration would lead 
to predictions necessarily absent from the simpler case. 

Let pi(e — j)'1d(e — j)'1 be the probability that the random 
fluctuation at s, has a value which lies between (¢ — j)'1 and (e — Jj)‘ 
- d(e — J)’; so that if 


erm (e = ya (1) 
V2 C1 
then, assuming a Gaussian distribution, 
1 
jy ote (2) 


x V/2n O71 
where o, is the standard deviation of the distribution function at $,. 
Similarly if 
BTS 
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Oy = (eI) (3) 
V2 02 
1 rach Fa ae (4) 
V2 02 


where o, is the standard deviation of the distribution function at 8, 
and p, is the ordinate of the same function. 

After the stimuli S, and S, have been simultaneously presented 
(Landahl, 1940), there develops at s, and s, the excitation factors & 
and e, respectively. These are functions of S; and S, respectively, and 
both are functions of time. Considering only the stationary state, so 
that, as t >> 1/a in equations (4) of the former development (A), 
the latter may be written, because of equations (1) of (A), 


scopes 

& pe 8 7. (5) 
ea TESS 

Eg urge ay ° (6) 


Superimposed on the amounts ¢, and «, at s, and s, respectively 
[Figure 1 (A)] are the random changes (e — 7)’, and (e — j)'2. 
The total amount at s, is then «, + (e — 7)’,, and that at s, is e + 
(e —7)'.. In order that S, produce its corresponding response, which 
we shall call “correct”, it is necessary and sufficient that the differ- 
ence between the total amount of the excitation factor at s, and the 
total amount at s, be greater than some quantity h (A, pp. 111, 115, 
and B, p. 165). Then the threshold, h’ , of fiber III of Figure 1 (A) 
is exceeded and a “correct’”’ response produced. For a correct response 
we have 


é:+ (e—Jj)4 &2 (e—J)'2>h. (7). 
Similarly for a wrong response we have 
éy + (e —J)'1— & — (e —j)'2< —h; (8) 


and for no response, which we shall identify with an “equality” 
(doubtful) response, we have 
h>eat (e—j)1— & — (e —J)'2 > —h. (9) 


Now the probability, f , that (e — 7)’, and (e — 7)’. each have 
specified values is given by f = p; p2 if the events are uncorrelated, 
but if not, then, designating the correlation coefficient by 7, we have 
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_ (a41-27 @142+a72) 
as 1 F (1-7?) » (10) 
220, 02V~1—?r 


Inequality (7) may be written 


(e —j)'1— (€ —J9)'2 >A—a te, (11) 
or defining 
(e—j)’ = V 2e0a= (e—j)1— (e—J)'2, (12) 
so that 
eh = 0; hy. — o> Oe ; (13) 


(11) may be written 
(e—j)' >— (a1 —& — A). (14) 


The probability, p(e« — 7)’ that (e — 7)’ has a particular value is 
given by replacing \/2 o2 a2 by V/ 2 01 a1 — V2 oa in equation (10) 
and integrating the resulting function of a,, f(a), with respect to 


\V/2 «, a, over all possible values of o,, that is, from —o to +o. 
Therefore, introducing (13) or its equivalent, a, = (0,0, — 00) /o2, 
into (10), and rewriting the exponent as a perfect square and a con- 
stant term, we obtain, on integrating. 


er a?/ (074-27 0102+ 072) 


a a SS Oe 
p (Vv e ) 220,02 V1—?” 

+00 

1, | 2 vorsriccoy - | = 
@ 3 1 VO71-27 0102+ 072 d ( V2 O1 a1) (15) 
-0 
e-772?/ (024-27 0102+ 072) 

SS eS 

V2 2(07, — 27 o; 02+ 072) 


or because of (12) 


1 a (=f) 
See Se 2(0%4-27O102+072) | (16) 
Van Vo7, — 27 01 02 O75 


Equation (16) is the equation of a Gaussian curve whose standard 
deviation is the square root of the expression in parentheses in the 
denominator of the exponent. We may for convenience define that 
value as o without violating the previous definition of of equation 
(12), introduced. to make the equations symmetric; or 


ple—j) 


o = Vo" — 27 0, 02+ o2 -<(17) 
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and equation (15) may also be written 


z Ci de (18) 


a“ V2ne0 


The probability of a correct response, P-, is given by integrat- 
ing p(e — j)’ over all values which satisfy the condition for a correct 
response given by (14) which is the same as (7). Then integrating 
p(e —j)' with respect to (e — j)’ from == (8; —€, —= Nh) AO oy we 
obtain 


Pa p(e—j)'d(e—j)’. (19) 
—(€1-E2-h) 


Similar expressions may be written for the probability of a wrong 
response, P,,, and the probability of an equality response Pz. Al- 
though each stimulus has a separate fluctuation associated with it, the 
mechanism predicts that only a resultant fluctuation will be observed, 
and it will be a function of the situation as a whole. 

Hence, if there are two independent or partially independent 
fluctuations instead of one, we obtain the same equations where a 
particular average of the two fluctuations, given by (17), is used. 
We have considered the separate variations as occurring at s, and s, 
of Figure 1 (A). Instead we may consider the random fluctuation as 
occurring at any point along the line of fibers traveled by the two re- 
spective groups of impulses (represented in the figure by fibers I and 
II). 

Consider the following case. Let the mechanism by symmetric 
so that we have o, = o2 and hence 


o=—o,V1—r. (20) 


Let a group of stimuli S, , S, --- be presented in pairs so that from the 
distribution of judgments we may obtain the standard deviation in 
terms of some function of the external stimulus, in particular, in 
terms of the logarithm of the stimulus (or what is equivalent, since 
S; © S,, in terms of the per cent change in the stimulus). The value o 
is measured in units of «. If k = a/AK , we obtain ke = ko’ , a partic- 
ular value of the standard deviation in terms of the per cent change in 
the stimulus. Now let the external conditions be changed such that 
we may assume that the parameters of the mechanism of Figure 1 
(A) are unchanged but that o, and hence o of (20) may have changed. 
We shall also assume that r of (20) is unchanged. Then let the same 
set of stimuli be presented in the same manner as before so that we 
may obtain ke = ko” , another experimental value. 
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Now let the stimuli be again presented in the same manner ex- 
cept that always one of the pair is under the condition for obtaining 
ke «’ and the other is under the condition for obtaining k o” . We then 
obtain a value ko . If the only change can be attributed to the changes 
in o, then o will be given by an expression of the same form as the 
right member of equation (17) with, however, o: replaced by 
o', = o' /V1—Tr and o replaced by o”1 = o"/V1—r [equation (20)], 
or 


¥ 
a ye ae ee a ae (21) 


V1l—r 


and hence 


ko = 


pe Sy (ko')? — 2r (ko’) (ko") +- (ko")?. (22) 
V1l—r 

The terms on the right of equation (22) are experimentally given ex- 
cept for 7. Since ko is an empirical value, 7 may be evaluated for 
any one set of values ko , ko’ and ko” . Choosing different experimental 
conditions, a new set of values ko , ko’ and ko” will be obtained and r 
may again be computed and so compared with the previous value of 
yr. However, it must be remembered that besides assuming 7 to be 
constant throughout, we assume that AK/a and h, of equations (5) 
and (6) are unchanged. Should one of the assumptions not be ten- 
able, we may change our purpose and use the same data to attempt 
measurement of the variation of one of the parameters. 

It was previously pointed out (B., p. 163) that, on the basis of 
the mechanism of Figure 1 (A) and on the assumptions made, a scale 
value was proportional to «,(S:) — é.(S2), both being proportional 
to the abscissa deviate corresponding to the ordinate of the cumulative 
frequency curve equal to the proportion of judgments SS, Se. Now 
o, and o. of equation (17) are measures of the degree of fluctuation 
associated with the stimuli S, and S2 respectively ; « measures the de- 
gree of the resultant fluctuation associated with the stimuli as they 
are compared, and is a measure of the distribution of the judgments 
of the pair. Equation (17) then gives the resultant fluctuation for a 
pair of stimuli in terms of that of each stimulus. If we measure 


-e,(S,) — &(S2) in terms of the standard deviation of the resultant 
distribution we may write 
€:(S;) 7 &2 (Sz) = X20; (23) 


but because of (17) 
€,(S1) sawn fg (53) — ree Vo7,— 27 v1 o2 + 07s (24) 
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where X,. is the deviate (Guilford, 1936, ch. VII) obtained when S; 
and S, are compared and o.. = o, the indices being introduced for 
definiteness. 

That equation (24) is similar to the equation of the Law of Com- 
parative Judgment (Thurstone, 1927) is not accidental. (Notice dif- 
ferences in the notation.) Statistically the problem here is essen- 
tially the same. The similarity is due to the assumption that the left 
terms are algebraically additive as well as to the assumption of nor- 
mality. It is a fundamental assumption in N. Rashevsky’s develop- 
ment (Rashevsky, 1938) that the excitation factor, as well as the in- 
hibition factor, is subject to algebraic addition. 

In equation (24) none of the terms except X,. is directly meas- 
urable. However, if both sides of equation (24) are multiplied by k 
as in (22), the left hand side becomes the per cent difference in the 
external stimulus magnitudes and on the right is the product of the 
deviate value X,. and the standard deviation in terms of the per 
cent change of stimulus value. All terms are then experimentally 
measurable. The equation only implies that the distribution of cor- 
rect or wrong judgments if plotted in terms of a per cent change in 
stimulus, or even better, in terms of the logarithm of the stimulus, 
will be accurately represented by the integral of a Gaussian distribu- 
tion curve (Guilford, 1936, p. 100). If such is not the case to a de- 
gree of accuracy demanded by experiment, then the assumption im- 
plied in equations (5) and (6) [equations (1) A] is not sufficiently 
accurate so that it is necessary to leave unspecified the ¢(S) relation 
(B, p. 163). 

Summary. It has been shown that the previous assumption 
(A, B) of a random fluctuation confined to one synapse instead of two 
was justified. Except in unusual situations the substitution for o 
from equation (17) generalizes the previous results to include a ran- 
domness associated with both stimuli. An experimental procedure is 
suggested in which it may be possible to obtain evidence for some of 
the assumptions made or to measure the variation of some external 
parameter with external conditions. 

II. Discrimination of stimuli having differences in more than 
one modality. Thus far we have always assumed that the stimuli to 
be compared differ in but one modality. Let us consider the case in 
which the stimuli differ in more than one respect. Let the two stimuli 
S and S, the latter being the standard, have the components S,, S2 
*++ Si, ++ Sp and S,, S. --- S,, ---S, in the modalities 1, 2, --- &, ++ 


p. We need consider only those modalities in which the two stimuli 
are different. 
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We consider two mechanisms, either of which may be used. First, 
let us assume that there is a discriminatory center represented by the 
mechanism of Figure 1 (A) but in which the fibers / and II of the 
figure represent respectively a number of fibers from the sense or- 
gans of various modalities. Let these fibers be at least p-fold so that 
at synapse s, [Figure 1 (A)] there corresponds to each S; a value 
e,(S;) for the excitation factor. The stimulus S then excites each of 
the fibers represented by I of Figure 1 (A) by the respective amounts 
S, --- Sp and the stimulus S excites the fibers represented by IJ of Fig- 
ure 1 (A) by amounts S, --- S,. Corresponding to each value S;, there 


is a value e.(S;) at synapse s., Figure 1 (A). If we define x by 
as D D (ate 
H Vil = ek (S:) —> fod (S;.) ’ (25) 
=1 k=1 


so that x is the difference between the total excitation factor at s, and 
the total at s. [Figure 1 (A)], then the equations (11) to (15) of B 
may be used if x of those equations is replaced by x of equation (25). 
Since some of the fibers of J and IJ may be inhibitory, certain e(S) 
should be replaced by 7(S) (where 7 is the inhibitory factor) or we 
may formally permit « to be negative in equation (25) to avoid com- 
plicating the notation, although « is essentially a positive quantity. 

Instead of the above described mechanism we may consider one 
slightly different. We may assume that within each modality there is 
a mechanism like that of Figure 1 (A). Then let the efferent ends of 
each of these mechanisms lead to a similar central mechanism. Let 
the fibers corresponding to a response to S lead to one synapse, $1, 
and those fibers corresponding to S lead to another synapse, s,. Then 
at s, the total excitation factor is proportional to the sum of all 
en (Sx) — @x(Sx) for which Sx > S, while at s, the total excitation fac- 
tor is proportional to the sum of all e,(S,)—ex(Sx) for which Si < Su. 
Hence the difference between the excitations at s, and s, is propor- 
tional to e.(S.) — éx(S;) summed over all k’s and thus is proportional 
to x of equation (25). 

The two mechanisms described give the same result. They would 
not necessarily do so if we allow for possible effects of the modalities 
upon one another. Such interplay is quite likely to occur. For the 
time we shall consider only the simpler situations described above. 

Since the modalities are not the same, the «(S) functions will 
vary with the modality, and hence the expressions & (Sx). These 
functions will increase monotonically with the S;, each S; being above 
some threshold value. Let us neglect the threshold values or assume 
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that they are nearly equal. The graphs of the functions then begin 
near the origin and for at least a small range increase linearly. Over 
this range «(S;) is proportional to S; or, say, equal to lx e(Sx), 
where « is independent of the modality and any constants are ab- 
sorbed in L;. On the other hand, if each function «,(S;) may be ap- 
proximated by a logarithmic relation over a rather wide range, or 


S 
ex (S;,) = Ly log i. ’ 


then, for approximately equal thresholds, we have e(Sx) proportion- 
al to log (S;/h) = e(S,), and « is again independent of k except 
through S. In either case and under even more general conditions, 
we may introduce the relation 


Ek (S;.) = Lye (Sx). (26) 


However, if the desired ranges in the S’s are not near the origins and 
some relations are linear and some logarithmic, etc., an expression 
similar to (41), but with a constant term added, would be a fair ap- 
proximation, or 


Ek (S:.) = lx é(Sx) + E0k « (27) 
Introducing either equation (26) or (27) into (25) we have 
= D D = 
L = = Ly, € (Sx) — > Ty, € (Sx) (28) 
or 
<= Dp 
x= = Ly Xx ; (29) 
where 
XL, = e(S,) — e (Sx) . (30) 


If we assume that the «(S) function is the same from individual to in- 
dividual so that any constant factors are absorbed in the L’s then 
e(S;), and hence x; , is a function of the stimulus only. The L’s vary 
from individual to individual and are the weights which each individ- 
ual accords to each modality. For example if two colored figures are 
being compared and it happens that one figure is accepted as being 
superior in form and the other is accepted as superior in coloring, the 
probability that an individual chooses one in preference to the other 
depends also on the relative weightings he gives to form and color in 
making the choice. 

Although the above results are sufficiently general to include com- 
plex stimuli, we have treated the stimuli as if they possessed only 
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physical attributes which could be measured. Some of the attributes 
which determine any judgment may not be physical. They may be en- 
tirely subjective and due only to a former association with the stimu- 
lus or a similar stimulus. Hence, some of the fibers leading to the cen- 
ter of discrimination arise directly from centers which have been con- 
ditioned to the stimuli in different degrees. If we give the proper set 
the subject may be able to inhibit all tendencies to respond on the ba- 
sis of the physical attributes so that (since we neglect non-influencing 
factors) all the factors, », [equation (29)] depend almost entirely 
upon the individual’s past experience. An illustration would be the 
ranking of nationality preferences by patriots of the various nations. 

Equation (29) has the form of the basic assumption in factor 
analysis (Thurstone, 1936) as the value x is proportional to the scale 
value of a certain comparison for a particular individual. If there are 
Q stimuli which are compared in all possible pairs (Guilford, 1936) 
there result Q(Q—1) /2 = N test situations. Let these N test situa- 
tions be presented to each of » individuals. By the scaling procedure 


(Guilford, 1936) we may obtain the value x;; for individual 7 on test 
situation 7. We then have the nN equations represented by 


= Dp 
Lj SET Tas (31) 
k= 


If the stimuli are chosen so as not to be dissimilar in too many re- 
spects, then, for any pair of stimuli, only a few of the scale values, 
Xj, , will be appreciably different from zero. Under these conditions 
we would be justified in applying the factor analysis methods devel- 
oped by L. L. Thurstone (1936). As it is only the correlations between 
the tests that are used in the analysis, one might substitute a very 
rough scale based on a simpler scaling method. 

In some cases the scale values of equation (29) may be deter- 
mined directly and, if so, would constitute a check on the factor pro- 
cedure and on the assumptions made. By modifying the factor proce- 
dure so that covariances are initially used, the actual magnitudes of 
the L’s and «’s could perhaps be derived in some common units. The 
process of correlating the results with information obtainable by 
other means should be simplified. It would also be of interest to no- 
tice the form of the distribution of the L’s. By careful selection of the 
items it might be possible to eliminate negative L’s — a negative L 
meaning an “unpleasant” attribute, being more strictly due to an in- 
hibitory process or “negative e”. 

It should be kept in mind that equation (29) would in general be 
an approximation to equation (25) and that certain simplifying as- 
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sumptions have been made. In particular it may be noted that no pos- 
sibility of intercorrelations among the k’s is allowed since in the group 
of fibers J and IJ [Figure 1 (A)] no interplay was considered. In gen- 
eral, one would rather expect some slight interrelation to be present 
in any single individual. But unless such an effect were fairly uni- 
form throughout the population it would not show up, but would con- 
tribute to the error space. The procedure suggested by equation (29) 
has been used as an empirical approach to the problem. The suggest- 
tion of the preceding paragraph may be of use in such an approach. 

Summary. By applying the theory to the comparison of complex 
stimuli, more general results are obtained. Under certain assumptions 
the results directly indicate the method of factor analysis. 

III. Applications to a multidimensional psychophysics. We have 
been considering discrimination between two stimuli which differed 
in more than one respect. From the very nature of the case we did 
not consider the possibility of the judgment being made from different 
points of view at different times. We might have assumed the more 
general situation in which each modality inhibited (or excited) every 
other, either as a function of the sum of the two components S; + S; 
or as a function of the difference, S; — S;,, or both. In this case, if 
either one modality was fairly prominent or if its difference was large, 
then the difference in magnitude in that modality might be the only 
factor affecting the judgment. 

In the simpler case having no such interrelation, we always have 
a difference between two sums of excitation factors which determine 
the judgment. In this case it would not be possible to give a spatial 
representation to a set of deviates corresponding to various stimulus 
pairs, the deviates being proportional to the differences between sums 
of «’s. However, a spatial represention in this case can be obtained if 
we separate the components of such a deviate into a factor which is a 
function of the individual and another factor which is a function of 
the stimulus pair. But this returns to the method of factor analysis 
described in the preceding section. Even in the more complex case of 
interrelation between modalities there would not be a result represent- 
able by coordinates of ordinary space, though, approximately, some 
curved surface may be used for representation of the data. 

A discussion of the possibility of representation of such data in 
spatial form has been given by Gale Young and A. S. Householder 
(1938). A rather elementary method for checking the possibility of 
spatial representation may be suggested here. Let stimuli be chosen 
varying in two modalities. Let them be chosen so that in many sub- 
groups they vary in one modality only. If the stimuli are represented 
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by Sj; then Si:, Sin, Sis, etc., have equal magnitudes in modality 7. 
With care in choosing the stimuli it would be possible to choose Say 
So; , Sa; , ete., so that these had equal magnitudes in modality 7. This 
is not necessary but simplifies the discussion. Since Si: , Sic} Diag Cle, 
vary in one modality only, then from the judgments involving these 
stimuli the corresponding deviates are obtained. Let us assume that 
any adjustment necessary is made so that the distance S,: to Siz plus 
the distance S,. to Si, equals the distance S,, to Sis. Choosing S,, as 
the origin and any direction to represent variation of j, the points 
Sar, Siz, Sas, ete., are determined. Similarly Sa. Saa Wo maCLCH are des 
termined along some other axis corresponding to?. Bute Sotponywosen 
etc., can also be determined. Thus all points S;; are determined (in 
fact highly over-determined). So far only variation of one modality 
at a time has been considered. 

We may now determine the distances involving differences in 
both modalities as S11 to Se2, or Si: to Ses, etc. The simple case con- 
sidered in the last section would require that the distance S,; to See 
would be equal to the sum of the distances S totSse ANGLS 1 tO. Searn lt, 
however, a spatial representation is possible, then the square of the 
distance S,, to So. will equal the square of the distance S,, to S12. plus 
the square of the distance Siz to So. minus the product of these two 
distances times the cosine of the angle between the axes. For each 
composite distance we can determine a value for the angle which will 
satisfy that one set of data. Thus, there are a large number of ways 
of obtaining the angle. If these values of the angle are consistent, 
then we may represent the data on a plane. If, however, the value of 
the angle varies significantly depending on the distances chosen, a 
plane surface cannot be used though some curved surface may be 
satisfactory. 

In order to obtain a psychological distance as discussed. above, 
it is necessary that there be uncertainty in the judgments. The range 
over which there is a sufficient degree of uncertainty is generally very 
small compared with the entire range. The ratio of the two would be 
approximately the reciprocal of the number of just noticeable differ- 
ences. In the above discussion of the possibility of a spatial represen- 
tation of the data, we have been considering only a small range of the 
space at a time. If it is necessary to use a curved surface to represent 
a small range of values, it would be of little use to attempt a spatial 
representation of the entire range. 

IV. Suggestions on the use of triads in psychophysics. The use 
of triads in presentation of stimuli for comparison judgments as sug- 
gested by M. W. Richardson (1938), has a number of advantages. 
However, in the manner of presentation used, a parameter is lost and 
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hence certain information is lost. In making a judgment from a triad, 
one is essentially making a comparison among three stimuli. The lat- 
ter may either be the three elements constituting the triad or the psy- 
chological distances between pairs, generally the latter. Now for such 
discrimination it was previously (A, p. 123) pointed out that one frac- 
tion of the correct judgments contain a false judgment. Hence, each 
of the respective proportions of judgments contain a certain fraction 
of false judgments and the situation is obscured. 

Consider the case in which the elements of the triad are A,B,C. 
Let the judgments be made on the basis of similarity so that a is the 
number of judgments that B and C are most similar or B is “closest” 
to C., ie. BC is the shortest psychological distance. Similarly let b be 
the number of judgments that, of the three elements, C and A are most 
closely related; and let c be the number of judgments that A and B 
are most closely related. We stated above that each of the number of 
judgments, as a, b, c, contain a fraction of false judgments. For 
instance, in deciding that AB is the shortest distance, a judgment con- 
tributing to the number c, the individual may have decided that BC 
was shorter than CA whereas the reverse may have been the case (as 
determined by the majority of individuals or by the majority of the 
same individual’s judgments). 

Now suppose that instead of permitting a choice from three re- 
sponses, the triad is presented so that the choice is limited to one of 
two responses. A particular triad would then have to be presented 
three times, as, AB compared with BC , BC with CA, and CA with 
AB. Let the triad be presented in such pairs a number of times. For 
the first pair, the comparison of the psychological distance AB with 
BC, let y, be the number of judgments AB > BC and let a, be* the 
number of judgments BC > AB. Similarly let a, be the number of 
judgments of BC > CA, and let £ be the number of judgments CA 
> BC in the second pair. Similarly for the third pair let 6, and V3 
be the numbers of judgments CA > AB and AB > CA respectively. 
Now it is the respective ratios which are constant, not the numbers, 
which naturally depend upon the number of presentations. Let us 
then multiply both numbers from the second pair, o and fb. by the 
fraction a,/a, and define 6, = f, o,/a.. Then the numbers a, and pf, 
may respectively be replaced by a, and f,. Similarly let us multiply 


both numbers from the third pair, 6; and y, by the fraction 61/Bs and 
define 


? o1 Be 
yee yee (32) 


Os 


* The a’s in this section are not to be confused with the a’s i i 
denote different things. SS ee 
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Then the numbers £3 and y; may be replaced by the quantities 6, and 
y',. Now y’, is a given number since a,, a, f., 63, and y; are num- 
bers experimentally obtained. The quantity y, is similarly an experi- 
mental value. Hence in general y, y',. There are then four indepen- 
dent values, as a,, 61, yi: and y’,. Other sets of four values are suf- 
ficient but there must be four. In the previously mentioned manner 
of presentation only three values would be obtained, a,b,c. Only if y, 
should be equal to y’, would it be possible to consider the similarity of 
the two sets of values. 

We next investigate conditions under which we should expect y’, 
of equation (32) to equal y,. For simplicity consider the differences 
AB , BC, and CA as being the stimuli being compared and let AB > 
BC > CA. Experimentally we obtain as above, yi, 01, 2, fo, Bs and 
ys. Then corresponding to the ratio y:/(y1 + a1) we obtain the prob- 
ability 


Passe = ee (33) 


that AB is judged greater than BC. Let Xy,0, be the deviate corre- 
sponding to Paprc . Then 


€y,a; — 91 Xy,2, C (34) 
Similarly, from o, and f2 
Oe 
= ——; 35 
Pacca Oo a Bo ( ) 
and Xz,g, is the deviate corresponding to Pacca , then 
Ea,8. — G2 X a8, . (36) 
And from #; and ys 
Bs 
Posas Bs =f Ys 


so that if X,,y, is the deviate corresponding to Posse , and Xg,., < 


0, then 
EBsys — G3 X Buys ? (38) 


As in equation (30) we may write 
Ey,04 OAR ——— ©BGL» 
Ea.B. — EBo — OA » (39) 


EBsYs = EoA — EAB- 
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The relations (39) hold provided we consider the simple case with no 
interrelation between the modalities, and if we assume no change in 
the parameters, in particular, no change in the relative weightings or 
“point of view”. From (39) it is evident that 


Ey,0; + Ea,B. Se EBsY3 — 0. (40) 


Hence, if in addition to the requirements just mentioned, we add the 
condition that the standard deviations are unchanged, or if 


Cigars OD eae OS (41) 
then from (34), (36), (38) and (40) 
Xy,0, rj Xap, + Xpiv5 =a (42) 


But when this condition holds, then, unless at least two of the dis- 
tances are equal, a very special case, it is impossible for y: to equal 
y',. That this is so can be verified by taking any numbers a > B 
> y and obtaining the deviates from the proportions y/(y + ©), 
a/(a-+ B) and 6/(6 + y). The sum will not be zero whereas by (42) 
the sum of deviates is zero. Hence the data, if they are represented 
by the mechanism, cannot be represented by the three values alone. 
We may also conclude that the proportions obtained according to the 
first discussed method, from which only three values a, b, and ¢ are 
obtained, cannot be used to determine corresponding deviates. 

It should be pointed out that, since e4g, esc and éc4 must be near- 
ly equal in a choice situation, if (40) and (41) hold, then (42) will 
be nearly correct. On the other hand, if (41) does not hold, or if there 
may have been some change in the parameters depending upon which 
“sides” of the triad are being judged, then (42) might not be even 
approximately true so that y, would not even be approximately equal 
to y',. Conversely, if y’, differs widely from y, in a case for which the 
distances AB , BC and CA are nearly equal, then the difference would 
be interpreted as due to a change in some parameters with a change 
in the pair of the triad to be judged. No attempt at this point will be 
made to give possible explanations on the basis of the mechanism of 
such a behaviour, should it occur. 

The author is indebted to Dr. A. S. Householder for reading the 
manuscript and for checking the algebra. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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: The generally accepted permeability theory of nerve conduction 
is presented in mathematical form. The resulting velocity formula 
is found to agree well with data on squid giant axon, but predicts 
velocities considerably too high in the case of Nitella. The depen- 
dence of velocity on fiber diameter is discussed for both medullated 
and non-medullated nerve, it being shown theoretically that velocity 
is proportional to the square root of diameter for non-medullated 
and to the diameter for medullated nerve. The equations relating 
the shape of the action spike to the observed permeability changes 
are given but are not solved. 


Conduction of the irritable response, while occurring in all types 
of living matter (Heilbrunn, 1938, chap. 38), has been most exten- 
sively studied, both theoretically and experimentally, in nerve. Three 
separate assumptions (verified in greater or lesser degree) have been 
fundamental in most discussion of the nerve conduction problem. 
These are: 

1. The hypothesis of local re-excitation by action currents, first 
put forward explicitly in 1872 by L. Hermann (1899). | 

2. The hypothesis of the nerve as a resistance-capacitance net- 
work; that is, that a nerve may be considered, as far as electrotonus 
is concerned, as a cable-like, linear electrical network. This point of 
view is sometimes attributed to J. L. Hoorweg (Cremer, 1909). 

3. The hypothesis of permeability change; that is, that excita- 
tion at a given point is characterized by an increase in membrane per- 
meability at that point. This is the basis of the J. Bernstein (1912) 
theory and has been extended by R. Lillie (1923). 

These three assumptions have come to be so widely accepted that 
diagrams of nerve conduction incorporating the first and third hy- 
potheses are given in the standard elementary physiology texts (e.g. 
A. J. Carlson and V. Johnson, 1937), while the second hypothesis has 
been used in almost all of the more exact treatments of conduction 
phenomena. In the present paper we shall give a quantitative account 
of nerve conduction which embodies these three fundamental hypothe- 
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ses; the related phenomena of excitation, and questions regarding 
the validity of these assumptions will not be considered here. 


I. Historical. To place the present theory in its proper perspec- 
tive we shall briefly recall some previous attempts to treat the prob- 
lem of nerve conduction quantitatively. One of the first of the mod- 
ern theories* was put forward by M. Cremer (1923), who obtained 
a formula for the velocity of propagation, but did not work out the 
details of the process. It appears that Cremer’s results are included 
under a later more detailed theory of N. Rashevsky (1938). K. Um- 
rath (1928), working with sensitive plants rather than nerves, ap- 
pears to have given the first quantitative account of a feasible con- 
duction theory involving a differential equation of propagation. His 
was a sort of “chemical fuse” theory, the excited region producing a 
substance which diffused ahead and ultimately brought about exci- 
tation at other points when its concentration there reached a certain 
critical value. 

Rashevsky (1931, 1933, 1935, 1988) considered the spread of ac- 
tion currents in a pure resistance network, used a formal excitation 
theory to describe the establishment of excitation by these currents, 
and represented the occurrence of excitation by the sudden insertion 
of a battery in the circuit at the point in question. By using different 
excitation equations he obtained slightly different conduction theories. 
He also considered myelinated nerve fibers, with excitation occurring 
only at the Ranvier nodes and conduction taking place by the distur- 
bance “jumping” from node to node. W. Rushton (1937) gave a simi- 
lar treatment, except that instead of using a formal excitation equa- 
tion he supposed excitation to occur when the membrane became suf- 
ficiently depolarized. His theory is equivalent to one of Rashevsky’s 
(Weinberg, 1940). 

A somewhat unsuccessful attempt to extend Rashevsky’s theory 
by including the membrane capacitance as an element in the Hoorweg 
network has been made by A. M. Weinberg, (1939). The problem can, 
however, be effectively handled by the methods used in the present 
article. 

The Rashevsky-Rushton theories introduce excitation formally 
as the sudden appearance of an E.M.F. in the circuit representing the 
nerve. In 1938 K. Cole and H. Curtis replaced this assumption by that 
of Bernstein; namely, that the essential event which marks the occur- 
rence of excitation is a decrease in the transverse resistance of the 
membrane. In their paper, however, these authors were not primarily 


* We are disregarding earlier theories of L. Hermann (1899 and B 
(1921, 1924) since they do not seem to conform to the Sere tipo, ‘of ews sae 
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interested in giving a complete theory of conduction, and so they did 
not attempt to derive a velocity formula. In the present work* we 
shall extend Cole and Curtis’s formulations with the purpose of obtain- 
ing a velocity formula and describing quantitatively as many details 
of the conduction process as possible. 


II. The Bernstein Model. Bernstein supposes the resting mem- 
brane to be polarized, i.e., to sustain a potential difference (corre- 
sponding to the resting or demarcation potential) like a charged con- 
denser. If the two sides of the membrane condenser be connected by 
a conductor there will be a passage of current and the charge will de- 
erease. If the conductor is removed the charge builds up to its original 
value, but this will not occur instantly. There is a time delay in the 
recharging process, as if the condenser is supplied from a battery 
through a series resistor. Specifically, Bernstein and Lillie picture 
the membrane as permeable to one ion of a solute but not to the other, 
so that a potential difference is set up at the membrane (Heilbrunn, 
p. 316; Erlanger and Gasser, p. 133). If the ions are removed from 
the surface it takes time for fresh ones to diffuse out from the interior, 
so that the resistance through which the battery supplies the con- 
denser may be essentially a diffusion resistance. The membrane is 
not completely impermeable to the other ion, so that the condenser is 
shunted by a leak resistance. Excitation is accompanied by a decrease 
in this shunting resistance. 

Picturing each infinitesimal element of the membrane in this way 
and linking the elements together by resistances corresponding to the 
longitudinal resistances of the interior and exterior regions of the 
nerve fiber, we have a network as shown in Figure la. Here 7; and ’ 


o) 
a> 


(b) 


FIGURE 1. 


are respectively the internal and external resistances per unit length ; 
C is the capacitance of the membrane per unit length; 1/p is the trans- 
verse conductance of the membrane per unit length; E is the voltage 


* The main results presented in this paper were obtained by F.O. in 1937 
and reported in abstract form (Offner, 1939) ; they were arrived at independently 
by G.Y. and A.W. in 1939. 
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of the battery which charges the membrane; and 1/F is the internal 
conductance of the battery per unit length. The total conductance per 
unit length, 1/s , between the interior and exterior of the nerve will 
be the sum of the membrane and battery conductances, and the sum 
of 7; and 7, will be denoted by M: 


s=Rp/(K-+ p), 
(1) 
M=r7,4+?r. 


The network of Figure 1a is the same as one discussed by Gerard 
(1930, p. 521), which is very similar to the one suggested by Cole and 
Curtis (1938), shown in Figure 1b. They considered that the total 
transverse resistance may be identical with the internal resistance of 
the battery. The two networks are equivalent (for the conduction 
problem) if H’ is equal to #s/R. The present picture however, ac- 
counts for the simultaneous change in both the measured membrane 
potential and permeability by assuming a change in the single para- 
meter p ; it is also slightly more general since it allows for a continual 
expenditure of energy (represented by the flow of current through & 
and p) to maintain a polarization of the membrane. This energy may 
be associated with a portion of the resting metabolic energy expendi- 
ture required to maintain the polarization. (Gerard, 1930). 

By a standard application of circuit theory (Cole and Curtis, 
1938) the differential equation satisfied by V , the potential difference 
between inside and outside of the core, is found to be 


“eV M oV M 
iq quite ih Teeny a eg 


ined (2) 
where x denotes distance along the nerve and t is time. 

According to the permeability hypothesis p, and therefore s, is 
supposed to decrease when excitation occurs.* Excitation will be 
assumed to take place when a given point is sufficiently depolarized; 
that is, when V at a point drops to some critical value V.. When exci- 
tation does occur, the membrane thus “breaks down” in the sense that 
its transverse resistance drops to a fraction of its resting value. A 
little later the membrane returns to its original semipermeability, s 
returning to its former value; the nerve then returns ultimately 
(though somewhat later than s) to its original state. The intimate 


nature of the breakdown-recovery cycle undergone by the membrane 


* More generally, all the quantities p, R, r;,7,,C, and E may change, but 


it is probably not too far wrong to think of them as all constant except p; Cole 
and Curtis (1939) find that C changes by not more than 2% during action in a 
giant axon of squid, and the change in M is probably very small also. 
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is not completely understood; but we may very roughly approximate 
the cycle if we suppose s to drop suddenly to a smaller value s,, remain 
there for a time z of the order of the absolute refractory period, and 
then rise suddenly to its original value s. This simple case will serve 
as an illustration, and we turn now to its consideration. 


III. U-Shaped Conductance Wave. Consider a disturbance mov- 
ing along a nerve with velocity v. An observer moving along with it 
would, on the supposition just made, see a contour for the total trans- 
verse resistance s as shown in Fig. 2. At the front end of the excited 


2A 


V eee an 

FIGURE 2. 
region the resistance drops from s to s» while at a distance vz behind 
this point it regains its original value.* If the propagation has reached 
a steady state with constant velocity, the voltage contour V will also 
be fixed and unvarying to the observer; i.e., in (2) we have a solution 
in the form 


V=V(«#—v21), (3) 


y = « — vt being the distance measured along the nerve from a point 
moving with the disturbance.** As depicted in Fig. 2, the disturbance 
partitions the nerve into three sections in each of which, under our 
present assumptions, (2) applies with constant coefficients. Equa- 
tions (2) and (3) give for V as a function of y, 

a?V M M 


dv 4 


*The experimental curve giving the time variation of s during activity 
(Cole and Curtis, 1939) suggests that a better approximation to the actual s 
contour would be obtained if s, is assumed to return to some value s,, 8» < § 
< s, instead of to its original value after the time 7. The details of the analysis 
are essentially the same as in the present case and so will be omitted for sim- 
licity. : 

% Umrath employed this same steady state solution while Rashevsky used 
an integral formulation and did not assume in advance that the velocity was 


constant. 
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which is the space differential equation for the steady V contour mov- 
ing with the disturbance in steady state conduction. The general so- 
lution of (4) is 


V =p E+ Acw + Bem (5) 


where w = a and w = f are roots (assumed distinct) of the quadratic 


equation 
wo? + vMCw — M/s = 0 (6) 


and A and B are integration constants. It may be noted here that a 
and £ have opposite signs; for definiteness, we suppose 8 > 0,a <0. 

Denoting the length of the depression in the s contour by 2 A and 
taking the origin for y at the center of this region, we have 


V,.=E'+A,e"4+ Bie (y> A) 
Ve SE, + Avewea Bo ele (—A ay <7) ee 
Vo EP Ager =P Brey = —— A) 


where LE’ = Es/R is the resting voltage at a point far ahead of the 
excited region and the zero subscript on a quantity denotes that s has 
been replaced by s, in the original quantity. EH’, is the voltage which 
would exist in the excited region far from the junction, if the break- 
down persisted indefinitely. 

To evaluate the integration constants, the following boundary 
conditions are used: 

a) V is finite at y = + o. 

b) By Kirchhoff’s laws, the voltage and current, and therefore V 

and dV/dy , are continuous at the junction points y = + A. 

From the first boundary condition, it follows immediately that B, = 
A; = 0, while the other boundary conditions will be satisfied if 


A, = (H'— E’,)J,/D 
A, = (E" — E',) J2/D 


(8) 
B, = (E' —E’,) J;/D 


B; = (E" — Ey) J,/D 
where 


J, = B (oo — Bo) e294 + B,(B — a) e-* 
—= (p Nast Bo) @ (2G0-2B0-a)A 


Jo = B(a — Bo) ce — a(B —B,) efa0-2804 
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J; = a(B — ao) eF4 — B(a— a) e264 (9) 
J, = (Bo — O%) ©8204 + a, (a — Bo) eF* 

— By (a — ay) @f24-2BorB)4 
D = (a— Bo) (B — a) — (a— a) (8 — Bo) era" P®. 
From (7), (8) and (9), the voltage at the front junction, V<, is 
ee ee (h, —— Eo) 
B (a9 — Bo) 7% + Bo (B — ao) — oo (6 — Bo) e?(2e Fo) 


6 ROS EE ee 
eae ae 
which for A large tends to 
¥.= fh — Bo CPLA AWE (11) 
Bo Ck 


It may be noted that the quantity V. is simply the voltage at the 
front junction between excited and unexcited regions. Bernstein sup- 
posed the membrane to break down when it had been sufficiently de- 
polarized — that is when V had reached some critical value V’.. If 
the breakdown is instantaneous, then V, and V’, are identical ; but, 
if as is probably the case, there is a time delay between excitation (or 
attainment of the critical voltage), and the breakdown of the mem- 
brane, then V, and V’. will be different. In the present theory only 
V., which can be determined directly from the point of inflection of 
the action spike is used, and so the calculation of the velocity is im- 
dependent of the excitation criterion used. In the discussion of the 
dependence of velocity on diameter the further assumption is made 
that V. (or more accurately E’/V.) for a given kind of nerve is inde- 
pendent of the diameter. 

The voltage at the point y = — A where the membrane resis- 
tance suddenly returns to its original value is given by an analagous 
expression which for A large tends to 


Qa 
V,=E'+ — (E'—E') . (12) 
The roots of (6) are of the form 
p=b—a 
(13) 
a=—b—a 


where b >a >0,a=4, and b, > b. In terms of these quantities 
(11) and (12) become 
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bo —a 
By hyseng tito Peay DE 14 
Vek, ee (E Eo) (14) 
and 
bo +a : 
et Da eee f 15 


Thus it is seen that, at least for A large, the voltage wave is not sym- 
metrical about the line y = 0 ; V is greater at the moment the mem- 
brane breaks down than at the moment it regains its original semi- 
permeability. ; 


IV. Velocity of Impulse. We have in equation (10) a relation 
from which to determine the steady state conduction velocity in terms 
of the nerve parameters. However, if we use the data on squid axon 
or Nitella, we find that if the excited region exceeds 0.6 cm., equation 
(11) is a very close approximation to the more exact equation (10), 
the error involved in neither case exceeding 1 part in 500. Since the 
measured crest time for 12 meter per sec. conduction is about 5 X 10-4 
sec., the excited region is certainly larger than 0.6 cm., and we may 
safely use equation (11); in other words, we may neglect the effect 
upon points ahead of the advancing wave due to the recovery of the 
membrane in its wake. This was done by Rashevsky and Rushton 
who simply considered the disturbance as leaving the nerve perma- 
nently altered behind it. This simplification tends to very slightly 
(less than 0.1%) overestimate the conduction velocity. 

From (14) we find 


bo—a 
Rah K (16) 
where 
K = (E' — V.)/(E'— E’,) . (17) 
From (13) and (6) it is seen that 
2a = vMC, 
b= Va?+ 1/L?, (18) 
1 
TE = M/s : 


After some algebraic operations involving repeated squarings, (16) 
gives 


[K/L? — (1— K)#/Le"]? 


4 Ka? CR en) ae (1 — K)?2/L,2 ’ 


(19) 
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or, from (18) 
bab Ao? § — A? So 
C\/ Ma dy 88)(A + do) AS ApS)” 


meneverg—eCh / Vie) — lis doi — (2 / Ve) * If Sh << 8». and-/A.and 
do are of the same order of magnitude, we may neglect /? 8s) compared 
with 4,28 , and A,»s) compared with 4s. Formula (20) then becomes 


Ast 
Up) 2 ee 
CAV M8y A. (A+ Ao) 
while if EL’) << V.sothati, wl, 


gen ened Soe saat 
"TCI Ms (LL) 


Equation (22) for the velocity may be compared with some ex- 
perimental data on squid giant axon given by Cole and Curtis (1939) 
and Cole and Hodgkin (1939). For an axon of diameter d = 500u, 
the following data are given: 


Capacity = 1 X 10-* fd/cem? so that C = ad X 10° 
15 1102 dems 

Membrane resistance during activity = 25 ohms 
cm? so that s, = 25/ad = 160 ohms cm. 

Specific resistance of sap = 34 ohms cm so that 7; 
= 1386/2 d? = 1.7 X 10* ohms/cm. 


External resistance in oil = 6.2 X 10 ohms/em. According to Cole 
and Hodgkin’s measurements, in sea water the characteristic length, 
L(L? = s/M) is 2.6 times the characteristic length in oil. A simple 
calculation on the assumption that only the external resistance is 
changed when the axon is immersed in sea water shows that 7, in sea 
water is negligible compared to 7; so that 


M =7; = 1.7 X 10* ohms/cm. 
The constant 4 can be approximated from the inflection point of 


the rising action potential spike. The permeability change occurs at 
the time when the action spike is inflected; the voltage at this point 


(20) 


(21) 


(22) 


* The steps leading to (20) are essentially the same if C is also allowed to 
change, say to C, ; the result in this case turns out to 
Ag? S— A? 8, 
9 == ——— (20 A) 
V Mss, XX, CO, (X FY Ao) (Ag 80”? + rs) 
where » = C,/C. If C= C,,” =1, and (20 A) reduces to (20). 
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may be taken as V.. (A slight correction [which is negligible if So ts 
sufficiently small] might be made since the recorded action spike max- 
imum differs from the negative of V by a quantity of the order of 
E’,). Using the squid action trace roughly corrected for amplifier 
response and electrode width (Cole and Curtis, 1939), 4 may be esti- 
mated to be about 1.1. With this value of 4, equation (22) gives 
vy = 25 meters/second, a value which agrees, within limits of experi- 
mental error, with Hodgkin’s (1939) value for velocity of conduction 
in squid axon immersed in sea water. 

In Nitella, the present theory predicts conduction velocities about 
100 times too high. This may be due to structural complications neg- 
lected herein (double membrane) as well as the relatively slow con- 
duction change. 

The effect on the velocity of immersion in oil can be calculated 
if the assumption is made that only M changes when the external me- 
dium is changed. For squid axon in oil M = 7.8 10* ohms/cm, 
while in sea water, M = 1.7 X 10 chms/em. The ratio between the 
two M values is 4.5 so that from (22), the velocities should be in the 


ratio 1 to 4.5 ; this is approximately the ratio found by Hodgkin 
for squid axon immersed in sea water and in air or oil. 

Equations (20-22) express the velocity of conduction as a func- 
tion of the electrical constants of the nerve. To determine how v va- 
ries with d, the axon diameter, certain assumptions must be made 
concerning the variation of the electrical parameters with d. The 
most plausible assumptions seem to be the following: 

a) The electrical parameters (membrane capacity, resistivity 
and membrane conductivity), when reduced to absolute units (e.g., 
membrane capacity per sq. cm.) are independent of the diameter. 

b) Aand 4, being dimensionless, are independent of d. From 
assumption (a) it follows that the electrical parameters which are 
measured per unit length of axon, vary thus:* 


Tf; oe. $,80,G alle a. 


For an axon in sea water (and probably in vivo), r. << 7; so that 
TUCO TF ce Ly 0 : 
Using these assumptions in equation (20), we find ~ 
v=kd (23) 


where k is a proportionality constant. The only data on non-myelin- 
ated fibers (Pumphrey and Young, 1939) indicates approximately 


_* The assumption that r; «1/d? is not accurate if the axon is so large that the 
action current density is not spread uniformly over a cross section. 
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such a dependence of velocity on diameter (v = kd” where n = 0.614 
mt .O21): 

The present model, although closest to a non-myelinated axon, 
may be applied to a medullated axon if the fundamental assumptions 
are made that action occurs only at the Ranvier nodes but that the 
nodes are so close together that a continuous network may be used to 
approximate the discontinuous one assumed to represent the actual 
nerve (Weinberg, 1939). From this point of view the resistance of 
the myelin is seen to be of no significance as long as it is very large; 
the myelin simply makes a discontinuous network out of a continu- 
ous one. The further assumptions are made: 

a) The capacitance of the myelin sheath is negligible compared 
with that of the membrane. This is probably not too far wrong in 
view of the relative thickness of the sheath compared with the mem- 
brane. 

b) The width of the Ranvier nodes is some function of d, say 
F(d), and the internodal lengths are some other function f(d). 

c) As in the non-myelinated case, the electrical parameters 
measured per unit dimension are independent of the diameter. 

Under these assumptions we find 


C,s « dX F(d)/f(d); M « 1/@? 
so that 


f(a) 
ies FD? (24) 
There is some evidence to show that f(d) = k'd™, where m © 1 (Ge- 
rard, 1931) ; unfortunately, as far as we know, no measurements to 
determine F'(d) have been made. If we make the simplest assump- 
tion, namely F(d) = constant, then we find 
Ud... (25) 

Evidence on the relation between conduction velocity and fiber 
diameter for myelinated nerve is conflicting. It is generally agreed, 
however, that v «< d” but values of m range from 1 to 2 (Erlanger 
and Gasser, 1937). The most recent work (Hursch, 1939) supports 
nha. 

It is unfortunate that in none of the experimental work have 
f(d) and F(d) been determined together with v — and it is clear 
that internodal length and width of the Ranvier nodes, as well as di- 
ameter explicitly, determine the explicit functional relation between 
velocity and diameter according to the present theory. Since the form 
of f(d) and F(d) may vary considerably from fiber to fiber, the pos- 
sibility exists that those authors who report m > 1 are working in a 
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range where internodal length and diameter are not strictly propor- 
tional, or nodal width is not independent of the diameter. 

The predictions of the present theory on variation of velocity 
with diameter are not peculiar to this particular theory. For example, 
the Rashevsky-Rushton velocity formula* 


gested. (26) 


where o is the dimensionless safety factor (Gerard, 1932; Tasaki, 
1939; Rushton, 1937) and L is the characteristic length, leads to the 
same general dependence on diameter as the present theory. In fact, 
probably any linear network theory of this type will lead to a similar 
v(d) function. 


V. The Membrane Cycle. The model discussed in this paper 
can be generalized if, in the original equation (4), the coefficients are 
any assigned functions of y, depending on how the circuit constants 
at a point of the nerve fiber change as an impulse passes by. Pre- 
sumably one would have to resort to numerical integration or to the 
use of machines such as the M.I.T. integraph in this case, since ana- 
lytical solution of (4) with coefficients arbitrary functions of y is not 
generally possible. Sufficient data for such a calculation are afforded 
by the work of Cole and his associates, but the labor involved is con- 
siderable unless a machine is used. 

However, we may invert the point of view and use the shape of 
the observed voltage contour or action spike to tell us something 
about the variation of the coefficients in (4). If we suppose only s to 
vary we can, assuming that the constant values of the remaining 


coefficients are known, determine s as a function of y by rewriting 
(4) as 


ia aay paciigety 
oy a anh as hy ae Oe 


where primes denote derivations with respect to y. 
From (29) we can calculate how s must vary in order to account for 
the observed action wave. Since derivatives of V are involved in (27), 
: *If our formula (22), which is simplified from the more complicated (20), 
is written in terms of the characteristic length, 
ye (22A) 
Cs,rXVIi+A 


and the two formulae are seen to be similar in structure, although the excited 


region constants appear in (22A). If (20) were used, the constants of both re- 
gions would appear. 
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this method requires rather accurate determination of the observed 
wave form; it is, so to speak, rather sensitive to variations in wave 
shape. The derivatives of V have been determined experimentally by 
K. Cole and H. Curtis (1939) and by B. Katz and O. Schmitt (1940). 

It can be readily seen from (27) that a symmetris s contour will 
not produce a symmetric V contour, for if V is symmetric the term 
in V is anti-symmetric while the remaining terms on the right are 
symmetric, so that 1/s is not symmetric. To produce a symmetric V 
wave, the value of 1/s would have to be greater on the leading side of 
the wave than at corresponding points on the trailing side. 

In comparing (27) with experiment, it must be remembered that 
the recorded action spike, Vic: , is related to V by 


ae eV (28) 


where V, is a constant, at least if monophasic potentials are recorded 
with one electrode on a distant portion of the nerve. Substituting 
this value into (27) 


E 
== (— 54-0 Ve + @) - (29) 


lipo 17 sents 1V een ke) (30) 


so that, aside from a scale factor, 1/s depends on the two parameters 
y and V,. The other constants can be determined directly by experi- 
ment. In the case of Nitella, v ~© 1 cm/sec so that vC << 1/M and 
(30) reduces to 


Aha (31) 


As was pointed out by Cole and Curtis, the theoretical 1/s and 
 V curves should rise more steeply than the experimental curves be- 
cause of the finite width of the recording electrodes. The exact correc- 
tion due to the finite electrode width is difficult to calculate, although 
an experimental determination of its order of magnitude has been 
made by these authors. Since 1/s is very sensitive to small changes 
_in V, an adequate comparison between theory and experiment is pre- 
cluded until accurate, corrected values of V and its derivatives are 
available. 
The authors wish to thank Dr. Carl Eckart for his assistance in 
the mathematical formulation and solution of the problem, and Pro- 


102 MATHEMATICAL BIOPHYSICS 


fessor R. W. Gerard for suggestions offered by him during several 
conversations. This work was aided in part by a grant from the Dr. 
Wallace C. and Clara A. Abbott Memorial Fund of the University of 


Chicago. 


LITERATURE 


Bernstein, J. 1912. Elektrobiologie, p. 95. Braunschweig: Vieweg & Sohn. 

Broemser, Ph. 1925. “Das Prinzip der Nervenleitung.” Zevt. F. Biol., 83, 355-398. 

Carlson, A. J. and Johnson, V. 1987. The Machinery of the Body. Chicago: Uni- 
versity of Chicago Press. 

Cole, K. and Curtis, H. 1939. “Electric impedance of Nitella during activity,” 
J. Gen. Physiol., 22, 37-64. 

Cole, K. and Curtis, H. 1938. “Electric impedance of squid giant axon during 
activity.” J. Gen. Physiol., 22, 649-670. 

Cole, K. and Hodgkin, A. 1939. “Membrane and protoplasm resistance in the 
squid giant axon.” J. Gen. Physiol., 22, 671-687. 

Cremer, M. 1923. “Die Ableitung der Cremerschen Formel fiir die Fortpflan- 
zungsgeschwindigkeit der Erregung im Nerven.” Quart. J. Exp. Biol. Supp., 
96-97. 

Cremer, M. 1909. Nagel’s Handbuch der Physiologie des Menschen. p. 910. 
Braunschweig: Vieweg & Sohn. 

Erlanger, J. and Gasser, H. 1937. Electrical Signs of Nervous Activity. Phila- 
delphia; University of Pennsylvania Press. 

Gelfan, S. and Gerard, R. W. 19380. “Studies of single muscle fibers.” Am. J. 
Psysiol., 95, 412-416. 

Gerard, R. W. 1931. “Nerve Conduction in relation to nerve structure.” Quart. 
Rev. of Borol., 6, 59-838. 

Gerard, R. W. 1930. “The response of nerve to oxygen lack.” Am. J. Physiol., 
92, 498-541. 

Gerard, R. W. 1932. “Nerve Metabolism.” Physiol. Rev., 12, 512. 

Heilbrunn, L. W. 1938. An Outline of General Physiology. Philadelphia: Saun- 
ders. 

Hermann, L. 1899. “Zur theorie der Erregungsleitung und der elektrischen 
Erregung.” Pfliig Arch., 75, 574-590. 

Hodgkin, A. L. 1939. “The relation between conduction velocity and the electrical 
resistance outside a nerve fiber.” J. Physiol. 94, 560-570. 

Hodgkin, A. L. 1938. “The subthreshold potentials in a crustacean nerve fiber.” 

' Proc. Roy. Soc. B., 126, 87-121. 

Hursch, J. B. 1939. “Conduction velocity and diameter of nerve fibers.” Am. J. 
Physiol., 127, 181-189. 

Katz, B. 1937. “Experimental evidence for a non-conducted response of nerve 
to subthreshold stimulation.” Proc. Roy. Soc. B., 124, 244-276. 

Katz, B. and Schmitt, O. H. 1940. “Electric interaction between two adjacent 
nerve fibers.” J. Physiol., 97, 471-488. 

Lillie, R. S. 1923. Protoplasmic Action and Nervous Action. Chicago: University 
of Chicago Press. 

Lillie, R. S. 1929. “Resemblance between the electromotor variations of rhythmi- 
cally reacting living and non-living systems.” J. Gen. Physiol., 18, 1-11. 

seedy F. 1939. “Circuit theory of nervous conduction.” Am. J. Physiol., 126, 


FRANKLIN OFFNER, ALVIN WEINBERG, AND GALE YOUNG 103 


Pumphrey, R. J. and Young, J. Z. 1938. “The Rates of Conduction of nerve 
fibers in cephalopods.” J. Exp. Biol., 12, 458-466. 

Rashevsky, N. 1931. “On the theory of nervous conduction.” J. Gen. Physiol., 
14, 517-528. 

Rashevsky, N. 1933. “Some physico-mathematical aspects of nerve conduction.” 
Physics, 4, 341-349. 

Rashevsky, N. 1935. “Some physico-mathematical aspects of nerve conduction. 
II. Physics, 6, 308-314. 

Rashevsky, N. 1938. Mathematical Biophysics: The University of Chicago Press. 

Rushton, W. 1937. “Initiation of the propagated disturbance.” Proc. Roy Soc. B., 
124, 210-243. 

Tasaki, I. 1939. “The electro-saltatory transmission of the nerve impulse and the 
effect of narcosis upon the nerve fiber.” Am. J. Physiol., 127, 211-227. 

Umrath, K. 1927. “iiber die Erregungsleitung bei sensitiven Pflanzen.” Planta, 
5, 274-324. 

Weinberg, A. 1939. “Nerve conduction with distributed capacitance.” J. Applied 
Physics, 10, 128-134. 

Weinberg, A. 1940. “The equivalence of the conduction theories of Rashevsky 
and Rushton.” Bull. Math. Biophys., 2, 61-64. 


